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1. Introduction 

Recently, the excitons in semiconductor have been the subject of the intensive theoretical 
and the experimental investigation. One of reasons is its technological impact on quantum 
electronics and photonics. Especially, people have studied some optical properties of the 
excitons in the confined quantum system , such as the quantum dot, the quantum wires, 
and the quantum well ^-^■ 

It is well known that if the emitted light field is in squeezing state, the fluctuation of 
one of its quadrature phase components will be reduced at the expense of an enlargement 
of the fluctuation in the conjugate partner. The quadrature phase component with reduced 
fluctuation could be applied to the low noise optical communications and the ultra-precise 
measurements. Many efforts have been dedicated to generate the squeezing state of the 
light fleld both theoretically and experimentally [^15] . 



In this paper, the emission of the high density excitons in quantum well under classical 
pump field is studied. We assume the quantum well is placed in the microcavity. The 
excitons interact only with a single mode of the cavity field. The reservoirs for both the 
excitons and cavity field have been taken into account. The stationary analysis of the 
system have been carried out by using the Fokker-Planck equation. The small fiuctua- 
tions approximation around stationary state is used to studying the quadrature squeezing 
properties of the output light field. 

2. Model and master equation 

We know when the density of the excitons become higher, the ideal bosonic model of the 
exciton is no longer adequate (In the case of GaAs Quantum well, the ideal bosonic model 
become inadequate when the density of the excitons exceeds 1.3 x lO^cm"^ |]16[). One 
way of dealing with these deviations is to introduce an effective interaction between the 
hypothetical ideal bosons. We consider the GaAs quantum well is placed in a microcavity, 
and the excitons in the quantum well are pumped by the classical light with the frequency 
uJb- Under rotating wave approximation, we have the Hamiltonian as following: 
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H = Ho + Hn^ + Hr^ + haT\ + ha+Ti + /ifor+ + /iS+Fs, (1) 



where: 



+ ihg{aH -ab+) + ih{Eb+e-"^'^ - E*be''^^^) , (2) 

ujc and ijjfe are frequencies of the cavity field and the excitons respectively, g is the coupling 
constant of the excitons and the cavity field, we assume it is a real number. The fifth term 
of the equation (2) is the pumping term. The pumping field is taken as a classical driving 
field. E is in proportion to its amplitude. a'^{a) is creation (annihilation) operator of 
the cavity field. b'^{b) is creation (annihilation) operator of the exciton. The exciton and 
photon modes are coupled respectively to their reservoirs with continuum modes which lead 
to dissipation. Hr-^ and Hr^ are the free reservoir Hamiltonians. The reservoir coupling 
may be considered as phonon scattering process in the case of exciton. But it is a simple 
single photon absorbing process for the cavity field. r^(r2) and r]^(ri) are the reservoir 
operators of the exciton and the cavity field with coupling constants included respectively. 
G is the coupling constant of the effective exciton-exciton interaction, it is a positive real 
number. 

The time evolution of the exciton and the cavity field system is described by the master 
equation of the reduced density operator. In Schrodinger picture , the master equation, 
which is based on the Born-Markoff approximation, is written as following: 

— = -t[Ho, p] + 7i(2apa+ - pa^a - a^ap) 

+ j2{2bpb+ -pbH-bHp), (3) 

where 71 and 72 are dissipative rates of the cavity field and the exciton. 

3. Classical equation and small fluctuation analysis 



By using two modes positive P representation |T7], we can convert the operator equation 
(3) into c-number Fokker-Planck equation: 
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where, a 
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OLx is eigenvalue for the coherent state of the hght field and 0^2 is 



eigenvalue for the coherent state of the exciton. We have transformed to the rotating 
frames of the cavity field and exciton by: 
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In terms of the positive P representation, ai(a2) and at (02) are not complex conjugate 
to each other. The above Fokker-Planck equation (4) has a positive semi-definite diffusion 
matrix leading a equivalent stochastic differential equation: 
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(6) 



i]i{t) are stochastic fluctuation forces witli zero means and satisfy tlie delta correlated 
function: 



(7) 



< Vi{t) > =0 

<vt{tHit')> =0 

<vtitHit')> = 6,,6it-t') 

The exactly solutions of the non-linear equation (6) are still difficult to be obtained. Here 
we will use the small ffuctuation approximation around the stationary state to reduce the 
non-linear equation (6) to a linear equation for the ffuctuations. A preliminary question is 
whether the stationary state of the system is stable against to small perturbation. To this 
end, we need neglect the ffuctuation fores and ffrst derive the stationary solution. Without 
the ffuctuation force, equation (6) becomes: 
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(9) 
(10) 

(11) 



Now ai{a2) and ai (at) ^^^ complex conjugate as can be seen from the equation them- 
selves. The stationary solutions a^ (ctf^) and a^ (ctt^) are determined by: 



-71 a? + g<^2 = 



-7i«r 



gaf 



-72a^ + E- i2Ga^^af - t/a? = 
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(12) 
(13) 
(14) 
(15) 



Then we consider the case that the system deviates somewhat from its stationary state, 
setting: 



(16) 



ai = ai + 6ai 

ai = ai^ + Saf 

Substitute the equation (16) into the equations (8-11) and only keep the fluctuations in 
the first order, we get the hnear motion equations for 5ai and 5af as following: 



— -Ja = —A5a, 
dt 



(17) 



Where, n2 = a2 '^^25 ^^^ 
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5^ 72 + iAGn2 i2Gd 
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Now we seek 
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g -i2Gat^^ 72 - z4Gn2 
for the solutions of the equation (17) of the form e^^. The eigenvalue A can be obtained 

from the equation 



\A-XI\ = 0. 



:u 



in which / is identity matrix. The eigenvalue equation deduced from the equation (18) is 



[(71 - A)(72 - A) + g'Y + UG'niij, - Xy = 



(19) 



The eigenvalues are : 



Ai,2 = ^{[71 + 72 + t2V3Gn2] ± \/(-7i + 72 + ^2^36-^2)2 - Ag^}, 



A3,4 = -{[71 + 72 - i2V3Gn2] ± 



-7i + 72-^2v^G'n2)2-V}, 



(20) 
(21) 



When all the four eigenvalues have the positive real parts, the stationary state {a'^,af^) 
is stable. That is, when any small deviation develops from q;° and af^, a^ and af will 
return to af and a^^. If the values of parameters (G, g 71, 72) are in the range to make 



the absolution value of the real part of v (— 7i + 72 ± i2^/?>Gn2Y — Ag"^ larger than 71 + 72, 



then a^{af^) are stable solutions. Now we consider the case of the stationary atate. we 
have linear equations for 6ai and 6a^: 
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We abbreviate this equation as following 
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—5d(t) = -Ad(t) + D^f](t), 
ot 



(22) 



(23) 



and will discuss the squeezing properties of the output field based on it in the next section. 



4. Quadrature component squeezing of the output field 



The fluctuation spectra of the cavity-field in the stationary state are the Fourier transform 
of the correlation functions < ai(t + t), aj{t) > in the positive P representation: 



Sij{uj) 



e-'"^^ <ai{t + T),aj{t) > dr, 



(24) 



where < ai(t + T),aj{t) >=< ai{t + T)aj{t) > — < aj(t + r) >< aj(t) >. Substituting 
the equation (16) into the equation (24), we have: 



/ + 00 
e"*"^" <6ai{t + T)6aj{t) > dr. 



(25) 



The fluctuation spectra will be obtained [0 as following: 



S{u) = {A + iujI)-^D{A^ - iujI)-\ 



(26) 



where / is the identity matrix and T stands for transpose. A and D are defined by 
equation(23). After somewhat tedious calculation, the fiuctuation spectrum of the cavity 
field will be obtained: 

5n(^) = -I ^ ^ {^G^nl{^l + uj^) + [g^ + (71 + iu){^^ + lu - lAGn^)] 

X [/ + (71 - ^w)(72 -i^- iAGn2)\], (27) 
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g(^2 2j^2 

Si2{uj) = S2i{uj) = =-^[9^11 + 7? 72 + u;'72], (2 



S22{uj) = i j--^{4:G'^nl{-ff + u^) + [g"^ + (71 + 20;) (72 + iu + i4:Gn2)] 

X [/ + (71 - «^)(72 - «w + ^4^712)]}, (29) 



where 



A= \[{^i + iuj){-f2 + ii^) + gT + l^G'nii-fi + iujy\\ (30) 

The creation and annihilation operators for output field cioutit) and a^ut(t) will be 
presented by two quadrature phase componentsX^ as following: 

acratit) = |[A°«*(t) +iX°"*(t)]e^(''-™) 



a+,(t) = l[A°"*(t) -2A°"*(t)]e-*(^-^*) 



(31) 

2,A°"*(t) -2A°"*(t)]e-^(^-^*) 

r2 is the frequency of the cavity field and 6* is a reference phase. The normally-ordered 
squeezing spectrum for each quadrature phase component of the output field may be ex- 
pressed by the fiuctuation spectra of the cavity fielde |T^ : 



: Sl^\u) : = 2j,[e-'''Sii{u) + e^'' 822(^0) + 821^ + ^12(0;)], (32) 

: 5!"*(a;) : = 271 [-6-2^^5^1(0;) - e^''822{uj) + 821(00) + 812(00)]. (33) 



If we choose the phase 



= 9,-^, (34) 



where 9i is the phase of the coherent state of the light field. The normal squeezing spectra 
of the quadrature componentsXi for the output field are: 



: ^r (^) ■■ = ^^[(/ + 7172)^ + uj' + 2Gnl{g% + 7^2) 

+ (7? + 72 - 2<?' - 12G2n2 + 2Gj2nl)uj' - UG'nhl] , (35) 

: S-^\u) : = ^^[-(/ + 7172)^ -u^' + 2Gnl{g'j, + 7^2) 

- (7? + 72 - 2<?' - 12G%1 - 2G^2nl)uj^ + l2G^nl^l] . (36) 

The most interesting part is the low frequency part. At a; = 0, the corresponding values 

are : Now we consider the case of w = , it corresponds to the resonance output. That is: 

: STiO) : = ^^[(/ + 7172)^ + 2Gn^(/7i + 7^2) - l2G\h!], (37) 



: ^°"(0) : = ^^[-(/ + 7172)^ + 2Gn^(/7i + 7^2) + UG^nhf], {31 



Ao 

Ao 
in which Aq denotes A(a; = 0). In the case 

{9' + 7172)' + 2Gnl{g'^i + 7^72) < UG^n^l (39) 

the A+ component will be squeezed at a; = 0. If if 

2Gn^(/7i + 7^2) + l2G'nhl < {g' + 7172)', (40) 

the X_ component will be squeezed at a; = 0. It deserves to be mentioned that when G = 0, 
namely the interaction between the excitons can be neglected, the fluctuation spectra of 
the output field are zero. This means that in this case the output field is a coherent light 
field as expected. 

5. Conclusion 

The quantum statistical properties of the field emitted by the high density exciton laser 
have been discussed. The conditions for the stability of the stationary state also have been 
given. The small fiuctuation approximation around the stationary state is made to deduce 
the fiuctuation spectra. We show that either the output qudrature phase components could 
be squeezed under certain conditions. 
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